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The structural and electronic properties of Si(335)-Au surface decorated with Pb atoms are studied 
by means of density-functional theory. The resulting structural model features Pb atoms bonded to 
neighboring Si and Au surface atoms, forming monoatomic chain located 0.2 nm above the surface. 
The presence of Pb chain leads to a strong rebonding of Si atoms at the step edge. The fact that 
Pb atoms occupy positions in the middle of terrace is consistent with STM data, and also confirmed 
by simulated STM images. The calculated band structure clearly shows one-dimensional metallic 
character. The calculated electronic bands remain in very good agreement with photoemission data. 
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I. INTRODUCTION 

One-dimensional (ID) atomic chains created on semi- 
conductor templates have recently attracted much at- 
tention due to new phenomena characteristic for a re- 
duced dimensionality The most spectacular exam- 
ples include a breakdown of Fermi liquid theory [2, 3] and 
Peierls metal-insulator transition [J]. A specially attrac- 
tive route to create such ID structures is a process of self- 
organization of atoms into very regular arrays of parallel 
metallic chains on stepped semiconducting or insulating 
substrates [5, 6]. In this case electrons near the Fermi 
energy are completely decoupled from the substrate due 
to band gap in electronic spectrum of the substrate. The 
chain atoms are bounded to the surface by low energy 
states which do not contribute to electronic properties of 
the system. 

The Si(335)-Au surface is one of the simplest ex- 
amples of high-index surfaces which stabilize the one- 
dimensional structures, and has been studied by number 
of techniques, including: reflection high energy electron 
diffraction (RHEED) [7], scanning tunneling microscopy 
(STM) [8i, angle-resolved photoemission spectroscopy 
(ARPES) d, fipl and first principles density functional 
theory (DFT) [ll^. In particular, the STM topography 
data show regular arrays of monoatomic chains separated 
by width of Si(335) terrace and a few nanometers long 
d, 0] . The ARPES spectra taken in the direction paral- 
lel to the steps show two highly dispersive bands crossing 
the Fermi energy {Ep)^ thus indicating one-dimensional 
metallic nature of the system [3]. The structural model 
of the Si(335)-Au surface has also been proposed [8] and 
confirmed later by the DFT calculations [ijj, which well 
describes all the available experimental data. 

Recently, other materials (Na, Pb, In) deposited on 
Au decorated Si (335) surface have also been studied 
[13, [H, [ii, [H, [ii, [13. Perhaps lead is the most inten- 
sively studied among them. Unlike Si(335)-Au surface, in 



* Electronic address: krawiec@ kft . umcs . lublin. pl| 



which gold substitutes Si atoms in the surface layer, the 
Si(335)-Au/Pb is a representative of new class of systems 
in which the deposited material forms one-dimensional 
structures adsorbed on-top of the surface. The lead de- 
posited on flat Si(lll)-Au(6x 6) surface forms well or- 
dered monoatomic layers [^], and is known to 
weakly interact with the Si substrate [2l|. Thus, in the 
case of Si(335)-Au/Pb structure, one can expect a more 
pronounced one-dimensional character of chains as com- 
pared to clean Si(335)-Au reconstruction. Indeed, the 
deposition of 0.28 ML of Pb on reconstructed Si(335) sur- 
face (0.28 ML gives exactly 1 Pb atom per Si(335) unit 
cell) also leads to one-dimensional objects on the sur- 
face. The STM topography shows a few nanometers long 
chains placed in between Si chains of original Si(335)- 
Au surface [14]. The photoemission spectra taken in the 
[110] direction (i.e. parallel to the steps), show a highly 
dispersive band crossing the Fermi energy and quite flat 
bands in the [112] direction (perpendicular to the steps), 
thus indicating clear one-dimensional character of the re- 
construction. A more detailed analysis of morphology of 
the surface and of the ARPES data is difficult or even 
impossible without appealing to DFT calculations. It is 
the purpose of the present work to propose a structural 
model of the Si(335)-Au/Pb surface and calculate corre- 
sponding band structure. 

Here I present a structural model of the Pb chains on 
Si(335)-Au surface, derived from total energy DFT cal- 
culations. It features single Pb atom per Si(335) unit 
cell placed near the Au chain. The Pb atoms are bonded 
to neighboring Si and Au atoms, forming monoatomic 
chain located ~ 0.2 nm above the surface. This picture 
is consistent with the STM topography data of Ref. , 
and also confirmed by simulated STM topography im- 
ages. On the other hand, the calculated band structure 
for the present model clearly shows one-dimensional char- 
acter, i.e. a strong dispersion in the direction parallel to 
the steps and their lack in the direction perpendicular to 
them, and remains in good agreement with the ARPES 
spectra of Ref. 0|. The rest of the paper is organized 
as follows. In Sec. [TTlthe details of calculations are pre- 
sented. The structural and electronic properties of the 
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clean Si(335)-Au surface are briefly discussed in Sec. Hill 
The structural model of Pb chains, simulated STM to- 
pography images and the band structure are presented 
and discussed in Sec. \IV\ \V\ Wl\ respectively. The in- 
fluence of step-edge buckling on electronic properties is 
discussed in Sec. IVIIi Finally, Sec. IVIIII contains some 
conclusions. 



II. DETAILS OF CALCULATIONS 

The calculations have been performed using stan- 
dard pseudopotential density functional theory and linear 
combination of numerical atomic orbit als as a basis set, 
as implemented in the SIESTA code 122|, 0,^, ^ . 
The local density approximation (LDA) to DFT [23] , and 
Troullier-Martins norm-conserving pseudopotentials [28t 
have been used. In the case of Pb and Au pseudopoten- 
tials, the semicore bd states were included. A double-^ 
polarizedfDZP) basis set was used for all the atomic 
species [23^, ^2^. The DZP utilizes two radial functions 
for each angular momentum and additional polarization 
shell. The radii of the orbit als for different species were 
following (in Bohrs): Si - 5.13 (35), 6.59 {3p) and 5.96 
{3d), Au - 4.39 (5d), 6.24 {6s) and 5.79 (6p), Pb - 3.56 
{bd), 4.68 (5/), 5.30 {6s) and 6.48 (6p), and H - 5.08 {Is) 
and 4.48 {2p). A Brillouin zone sampling of 24 nonequiv- 
alent k points, and a real-space grid equivalent to a plane- 
wave cutoff 100 Ry (up to 82 k points and 300 Ry in the 
convergence tests) have been employed. This guarantees 
the convergence of the total energy within ~ 2 meV per 
atom in the supercell. 

The Si(335)-Au/Pb system has been modeled by four 
silicon double layers and a vacuum region of 18 A. Ah the 
atomic positions were relaxed except the bottom layer. 
The Si atoms in the bottom layer were flxed at their bulk 
ideal positions and saturated with hydrogen. To avoid 
artiflcial stresses, the lattice constant of Si was flxed at 
the calculated value, 5.41 A. The atomic positions were 
relaxed until the maximum force in any direction was less 
than 0.01 eV/A. 

III. CLEAN SI (335)-AU SURFACE 

The deposition of 0.28 ML of gold on Si(335) surface 
forms well ordered arrays of chain structure. The surface 
consists of (111) terraces which have a width 3| x a[xi2] 
(1.26 nm) [7]. Each terrace contains a single row of gold 
atoms running parallel to the step edge, i.e. in the [110] 
direction. The gold chain is formed by substitution of Si 
atoms in the middle of terrace. The step edge Si atoms 
form a 'honeycomb' substructure [29, 30, 31], which is 
a common feature of all the Au-induced vicinal Si sur- 
faces m. The structural model of the Si(335)-Au surface 
is shown in Fig. [TJ where the Si surface atoms (Sii- 
Sie) are labeled by numbers 1-6, and the gold atom by 
Au. Originally this model has been proposed as a sim- 




FIG. 1: (Color online) Structural model of Si(335)-Au surface. 
Top panel shows side view of the structure, and bottom panel 
shows top view with marked surface unit cell. Labels 1-6 
stand for silicon surface atoms (Sii-Sie), while Au denotes 
gold atom. The dashed lines in bottom panel indicate step 
edges. 

pie truncation of Si(557)-Au reconstruction [s], and later 
conflrmed by DFT calculations 111). 

The STM topography data [8, 9] show a single chain 
per terrace, which is associated with the step edge Si 
atoms rather than with the Au chain. The gold substi- 
tutes some of top layer Si atoms in the middle of terrace, 
and is not visible to STM. This has also been conflrmed 
by DFT calculations [11]. See also simulated STM topog- 
raphy images presented in Fig. [3l The ARPES spectra 
show two highly dispersive bands crossing the Fermi en- 
ergy in the direction parallel to the steps, and quite flat 
bands in the direction perpendicular to them [8], clearly 
indicating one-dimensional metallic nature of the system. 
The DFT calculations indicate that one of the bands 
crossing Ej? is associated with the step edge Si atoms, 
having unsaturated dangling bonds, while the other one 
originates from hybridization of the Au and neighboring 
Si atoms in the middle of terrace in the surface layer [JjJ . 



IV. STRUCTURAL MODEL OF PB CHAINS 
ON SI(335)-AU 

One-dimensional structures of Pb on Si(335)-Au re- 
construction one obtains assuming single Pb atom per 
Si(335)-Au unit cell, which corresponds to the experi- 
mental Pb coverage - 0.28 ML [14]. The total energy 
calculations show that Pb atoms prefer to adsorb on the 
surface. The substitution of Pb atoms into top Si layer 
is energetically less favorable, as the surface energy is 
by 0.3-0.6 eV (per unit cell) higher in comparison to 
clean Si(335)-Au surface and Pb atom in the bulk fee 
structure. Similarly, the substitution into the second Si 
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TABLE I: The relative surface energies of most stable struc- 
tural models of Si(335)-Au/Pb structure. The energies are 
referred to clean Si(335)-Au surface and Pb atom in the bulk 
fee structure. 



model 


position of Pb 


surface energy (eV) 


1 


Sis-Au-Au 


-0.41 


2 


Sii (subst) 


-0.31 


3 


Sis-Sie-Sii 


-0.25 


4 


Au-Sii-Sii 


-0.06 


5 


Si4-Au 


-0.02 



layer is not preferred. In this case the energy cost is 1 
eV. The exception is the substitution of Pb at the step 
edge, where the energy gain is 0.31 eV. More than 40 
structural models have been investigated, and only five 
of them lead to stable structures. Corresponding surface 
energies, with respect to clean Si(335)-Au surface and 
bulk Pb atom, are shown in Table [H The differences in 
energy are rather small, however, as it will be argued 
later, the model with the lowest energy, i.e. Sia-Au-Au 
model, is the best candidate for a true model of Pb chains 
on Si(335)-Au surface. Note the nomenclature used here 
reflects the bonding of lead with corresponding surface 
atoms. For labeling see Fig. [2l The structural model 
with the lowest surface energy features Pb atoms located 
2 A above the surface near the gold chain, and is shown 
in Fig. [21 As one can read off from the figure, lead is 




FIG. 2: (Color online) Structural model of Si(335)-Au/Pb 
surface. Top panel - side view, and bottom panel - top view 
of the model. Labels 1-6 stand for silicon surface atoms (Sii- 
Sie), while Au denotes gold and Pb - lead atom. Surface unit 
cell is also indicated in the bottom panel. 



bounded to two Au atoms and one Sis atom, but not to 
Si4 atom (Fig. [2]). The Pb-Sia bond length is equal 2.75 
A, and Pb-Au - 2.92 A. Both bonds are shorter than 
the Pb bonds in the bulk fee structure (3.46 A). The 
presence of Pb atoms on-top of the surface also modi- 
fies positions of underneath Au and Si atoms, leading to 
change in their bond lengths. Thus the Au-Si4 bond is 
equal 2.52 A, and Au-Sis - 2.47 A, to compare with 2.43 
A and 2.38 A, respectively, in the clean Si(335)-Au sur- 
face. What is more important, presence of Pb atoms also 
leads to a strong rebonding at the step edge. In the clean 
Si(335)-Au surface, the Si atoms near the step edge form 
a sort of honeycomb chain (see Fig. [T]), which is common 
feature of all the vicinal Si surfaces [8]. Main feature 
of this substructure is a true double bond between Si2 
and Sis atoms, which is responsible for the stability of 
the honeycomb chain [2^. In the presence of Pb atoms 
on Si(335)-Au surface, there is no longer double bonds 
between Si2 and Sis atoms. The atoms are just single 
bonded, as in common silicon structures. The Sis atom 
is now bonded with the lead, and Si2 with underneath Si 
atom in the second layer (compare Fig. [T]and Fig. [2]). 

At this point I would like to comment on main struc- 
tural features of the other models, listed in Table [U which 
have slightly higher relative surface energies. The next 
'best' structural model features the step edge Si atoms 
substituted by lead (model 2 in Table [Tj). Model 3, hav- 
ing relative surface energy -0.25 eV, accounts for Pb atom 
bonded to Sis and Sie atoms within the same terrace and 
Sii atom at the step edge of the neighboring terrace (see 
Fig. [2] for labeling). In model 4, with the energy -0.06 
eV, Pb is located between two Sii atoms at the step edge 
and one gold atom on neighboring terrace, while in model 
5 the Pb takes bridge position between Si4 and Au atoms 
on the same terrace. However those models can be ruled 
out, according to arguments given below. 

The present model, i.e. the model with the lowest rel- 
ative surface energy, is a good candidate for structural 
model of Si(335)-Au/Pb reconstruction. There are few 
strong arguments supporting this model. First one is 
that the model has the lowest energy. The second one is 
that the model explains the STM data of Ref. [14]. In 
STM images of Si(335)-Au/Pb surface, the Pb chains are 
located in between the chain structure of original Si(335)- 
Au reconstruction (see Fig. 1 of Ref. According 
to the STM data, the Pb chains should be located in 
the middle of terraces, above the surface. This is exactly 
what the present model shows, and STM simulations, 
which will be presented in the next section, confirm this. 
The other models, listed in Table [J remain in disagree- 
ment with the STM data. Model 1 leads to a single chain 
in STM topography, as the step edge Si atoms, which 
formed the monoatomic chains in clean Si(335)-Au sur- 
face are now replaced by the lead atoms. Similarly, mod- 
els 2, 3 and 4 also lead to single chain within terrace, as 
the Pb atoms saturate the step edge Si dangling bonds. 
In fact, the STM topography of model 4 shows a zig-zag 
chain, associated with Pb atoms and less visible step edge 
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Si atoms. Only the last model, i.e. model 5, remains in 
reasonably good agreement with STM topography data. 
The next argument concerns the band structure, as the 
only present model reproduces the photoemission spec- 
tra of Ref. [14] very well. The other models, which have 
slightly higher relative surface energies disagree with the 
ARPES data, in particular, they do not give correct band 
structure near the Fermi energy. This will be further 
discussed in Sec. |Vll The next argument, albeit intu- 
itive, is based on the analogy with the growth of Pb on 
fiat Si(lll) surface. It is well known, that lead grows 
in very regular fashion, i.e. layer by layer, on Si(lll)- 
Au(6x6) surface from very beginning [H, [l^, [iO], con- 
trary to the growth on Si(lll) with 7x7 reconstruction, 
where it forms amorphous wetting layer [18*, 's^, 'SS] . This 
indicates that Pb prefers to bond with Au rather than 
with Si atoms. The present model of Si(335)-Au/Pb sur- 
face also reflects this fact, as the Pb atom is bonded with 
one Si and two Au atoms. The last argument comes from 
the conditions at which the experiment of Ref. 0] has 
been performed. Namely, the temperature of deposition 
of Pb on Si(335)-Au surface was 260 K. This temperature 
is far to low to substitute the Si atom by the lead, and 
the model with Sii atoms replaced by Pb can be ruled 
out. All the above arguments show that the model shown 
in Fig. [2] is very good candidate to be a true model of 
the Pb chains on Si(335)-Au surface. 



V. STM SIMULATIONS 

The STM topography data of clean Si(335)-Au sur- 
face shows one-dimensional structures which are inter- 
preted as the step edge Si atoms jl, [ll| . The deposition 
of 0.28 ML of Pb leads to monoatomic Pb chains located 
in the middle of the Si(335)-Au terraces, i.e. between 
the Si chains of original Si(335)-Au surface. The struc- 
tural model discussed in previous section supports this 
scenario. To further check the validity of the structural 
model, I have performed STM simulations within the 
Tersoff-Hamann approach [34J. The results of constant 
current topography for different bias voltages are shown 
in Fig. [3l Top panels represent simulated STM topogra- 
phy of 3 X 3 nm^ of the same area of the Si(335)-Au/Pb 
surface for sample bias U = — 1 V (a) and U = —0.1 V 
(b). For comparison, corresponding images of the same 
area of clean Si(335)-Au surface are shown in the bot- 
tom panels. As it was discussed previously (see Sec. IIIip . 
the STM topography of clean Si(335)-Au surface features 
monoatomic chains, which are associated with the step 
edge Si atoms (see Fig. [3b) and d)). The Pb atoms de- 
posited on this surface form monoatomic chains, located 
in the middle of terraces, i.e. between the chains of clean 
Si(335)-Au surface. This is evident, if one compares pan- 
els a) and c) or b) and d) of Fig. [3l 

As one can read off from Fig. [3l the Pb atoms are 
more pronounced than the step edge Si atoms, especially 
at high sample bias. Such a behavior can be explained 




FIG. 3: (Color online) STM simulations of 3 x 3 nm^ area of 
Si(335)-Au/Pb surface (top panels) and Si(335)-Au (bottom 
panels) for sample bias U = — 1 V (a and c) and U = —0.1 V 
(b and d). All the images show the same area of the surface. 

by the combination of structural end electronic effects. 
The Pb chain sticks out (by 2 A ) above the surface, and 
this mainly contribute to the discussed effect. On the 
other hand, the electronic properties also play significant 
role. The calculated projected density of states (PDOS) 
of the Pb atoms and the step edge Si atoms is shown in 
Fig. m The Pb PDOS features larger values at ener- 
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FIG. 4: (Color online) Projected density of states of the Pb 
(solid line) and the step edge Si atom (dashed line). 

gies — 0.7 eV and +1.5 eV (see solid line in FiglJj), while 
PDOS of the step edge Si atoms is larger near the Fermi 
energy (dashed line). In first approximation, the STM 
current is proportional to integrated density of states be- 
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tween the Fermi energy and applied bias voltage (elJ). 
This contributes to the fact that the Pb chain is more 
pronounced at higher sample bias. At very low sample 
bias {U = —0.1) V, the step edge Si atoms have larger 
values of PDOS, and as a result both chains have com- 
parable topography amplitude (see panel b) of Fig. [3]). 
However, the Pb chain is still more pronounced than the 
Si one due to the fact that it sticks out above the surface. 

The STM topography along the Pb chain shows very 
small changes of the amplitude at low voltages (see panel 
b) of Fig. [3]), indicating strong overlapping of the 6p 
states. This is also reflected in the band structure, where 
6p band of lead, crossing the Fermi energy, is very dis- 
persive in the direction along the chain (see Fig. [4]). On 
the other hand, there is no such an effect in the STM to- 
pography of the step edge Si atoms due to more localized 
character of 3p states of silicon. As a result the step edge 
Si band is rather flat, i.e. is less dispersive than 6p band 
of lead (see Fig. [4|). All this shows that the Pb chain 
has more metallic character than the step edge Si chain, 
which makes the lead a good candidate to look for ex- 
otic phenomena characteristic for the systems of reduced 
dimensionality. 

VI. BAND STRUCTURE 

The calculated band structure for present structural 
model, along the high symmetry lines of two-dimensional 
Brillouin zone (BZ), is shown in Fig. O The direction in 
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FIG. 5: (Color online) Calculated band structure along high 
symmetry lines in two-dimensional Brillouin zone (shown on 
the right). The [110] direction is parallel to the steps, while 
[112] - perpendicular to them. The bands associated with Pb 
chain is marked with filled triangles, the step edge Si atom 
band - open circles, while the bands originating from the Au 
chain hybridizing with Si4 and Sis atoms are marked with 
squares. The Fermi level is taken as the zero of energy. 

the 2D Brillouin zone, defined by the points F, K and M', 
is parallel to the steps of the Si(335)-Au surface, while the 
F-M is perpendicular to them (see right panel of Fig. [5]). 
A few important surface bands are marked with symbols 
in Fig. [5l The band marked with open circles pinning 
the Fermi energy, originates from unsaturated bonds of 



the Si atoms at the step edge. The band 0.5-1.5 eV be- 
low the Fermi energy (open squares), as well as more 
dispersive band crossing the Ef (filled squares) are as- 
sociated with the Au chain. To be more precise, those 
bands originate from the hybridization of the Au chain 
with neighboring Si4 and Sis atoms (see Fig. [2]). All 
the bands discussed above have also been identified in 
the clean Si(335)-Au surface [11] and similar bands have 
been found in the Si(557)-Au reconstruction In 
fact, the band marked with filled squares has also Pb 
character. This band reflects the hybridization of gold 
with Pb and Sis atoms. Finally, the band shown as filled 
triangles is associated with 6p states of lead. All the 
above bands do not have (or have very weak) dispersion 
in the direction perpendicular to the steps (F-M), indi- 
cating one-dimensional character of the structure. 

Since Pb features a strong spin-orbit (SO) interaction 
HH, [33 , and the low-dimensionality of the system can in- 
crease it [38], it is worthwhile to comment on this effect 
in the Si(335)-Au/Pb system. Although the spin-orbit 
interaction was not included in present calculations, one 
can draw some conclusions appealing to the Si(557)-Au 
reconstruction. The measured band structure of Si(557)- 
Au surface shows two Au induced proximal bands cross- 
ing the Ef near the K point (of 2x1) zone with ~ 300 
me V splitting |l3, En] . An explanation of the split- 

ting was given in terms of SO interaction [42]. Since both 
Si(335)-Au and Si(557)-Au surfaces belong to the same 
family of vicinal surfaces (Si(335)-Au reconstruction may 
be considered as a simple truncation of Si(557)-Au sur- 
face), one can expect SO interaction to play similar role 
here. In fact. Grain et al. [8|] observed small splitting 
of the band crossing the Ef- This was interpreted as 
two different bands, one coming from Au-Si hybridiza- 
tion and the other one from the step edge Si atoms [ll[. 
It is also possible that the splitting has its origin in the 
SO interaction, similar like in Si(557)-Au surface. On 
the other hand, if Pb is deposited onto Si(335)-Au sur- 
face, one would expect SO effect to be more pronounced. 
However, the ARPES data [14] shows a single metallic 
band and thus no evidence of SO splitting. So one can 
conclude that if SO effect is really present in this sys- 
tem, certainly the splitting of the band is smaller than 
the energy resolution of the ARPES apparatus, which is 
50 meV in this case. This is rather an intriguing result, 
however it can be true, as recent spin-polarized photoe- 
mission measurements report values of SO splitting as 
small as 15 meV in Pb deposited on Si(lll) surface [4^. 

A comparison of the calculated band structure with the 
photoemission spectra of Ref. [14] is shown in Fig. [6l As 
one can see, the present DFT calculations agree very well 
with the experimental data. In particular, the band near 
the Fermi energy, which has high photoemission inten- 
sity, is the 6p band of lead. One can also identify other 
bands, marked with symbols in Fig. O except the band 
coming from the step edge Si atoms. The reason that this 
band is not observed in ARPES can be twofold. First, 
the band is to close to the Fermi energy which makes 
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FIG. 6: (Color online) Comparison of the measured ARPES 
intensity along F-K-M' line of the 2D Brillouin zone (Ref. 
p^ ) and the calculated band structure. The high intensity is 
shown light. 

it very difficult to verify experimentally. Second, there 
may be some mechanism which leads to the saturation of 
the Si bonds, accounting for various imperfections or im- 
purities. It is more likely that the band associated with 
the unsaturated step edge bonds (open circles in Fig. [5]) 
is eliminated in the real surface by some reconstruction 
and passivation. The main effects of the reconstruction 
can be achieved by saturating the dangling bonds with 
hydrogen [35] . This process removes extra band from the 
gap region (open circles in Fig. [5|), and is energetically 
very favorable, as the energy gain is ~ 1.30 eV with re- 
spect to the unsaturated surface and H2 molecule. The 
comparison of the calculated band structure with satu- 
rated Sii bonds and the ARPES spectra of Ref. [14] is 
shown in Fig. [71 Clearly, the band structure is in bet- 
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FIG. 7: (Color online) Comparison of the measured ARPES 
intensity (Ref. [ij) and the calculated band structure for the 
model with saturated step edge Si bonds. 

ter agreement with the ARPES experiment now. This 
suggests that the step edge Si bonds are really saturated 
by some surface reconstruction, however this cannot be 
accounted for without multiplying the unit cell of the 
Si(335)-Au surface. 



VII. STEP-EDGE BUCKLING 

The Si(557)-Au reconstruction is known to undergo a 
buckling of the step edge at low temperatures [36, 42, 44]. 
The step edge Si atoms alternate between up and down 
positions, with distortion in z direction = 0.65 A. In 



general, all the step edge atoms occupy equivalent posi- 
tions and feature half- filled dangling bonds, which give 
rise to a flat band near the Fermi energy. This is re- 
ferred as high- temperature phase. However this situa- 
tion is energetically unfavorable at low temperature and 
the system tends to lower its energy by the buckling of 
the step edge. The energy gain is 130.9 meV. The dan- 
gling bonds of the up-edge Si atoms became fully occu- 
pied, while those of the down-edge are empty. As a result 
the band structure features the fully occupied and empty 
step-edge electron bands. The buckling of the step edge 
influences the properties of the band structure near the 
Fermi energy through changes in the Au-Si-Au bond an- 
gles. Those changes are responsible for the band gap that 
opens in dispersive Au-Si band crossing the Ei?, and thus 
for the metal-insulator transition at low temperature ^43 |. 

Doubling the unit cell of Si(335)-Au/Pb system in di- 
rection [110] one can account for similar effect associated 
with the buckling of the step edge. In the present system 
the distortion is slightly smaller than in Si(557)-Au 
system, and equal to 0.57 A. This smaller value of the 
Az can be related to the fact that Si(557)-Au supercell 
cell consists of two non-equivalent x 1 unit cells with ex- 
tra Si adatom in one of them, which in turn can increase 
Az. The x2 Si(335)-Au/Pb supercell is constructed from 
two equivalent xl cells. The buckling of the step edge 
lowers the total energy of the system by 117.5 meV, to 
be compared with 130.9 meV in the case of Si(557)-Au 
surface, and leads to a splitting of the step edge Si band, 
originally pinning the Fermi energy. As a result one ob- 
serves two bands, a fully occupied flat band (with a band- 
width W ^ 0.09 eV) at energy 0.29 eV below E^, and 
empty more dispersive band {W ~ 0.38 eV) at energy 
0.25 eV above the Ep. Those bands are associated with 
up-edge and down-edge Si atoms, respectively. Similar 
bands have also been observed in the Si(557)-Au surface 
[4^ ]. The step edge buckling is also reflected in different 
bond angles between Au and Sis atoms (see Fig. [2]). In 
situation when all the step edge Si atoms occupy equiv- 
alent positions {Az = 0), these angles are equal 101.64°, 
to be compared with 111.6° and 103.7° for Si(557)-Au 
surface. Again, different values of the Si(557)-Au angles 
come from non-equivalent x 1 unit cells, and the buckling 
can change them by ±10° [3]. On the other hand, the 
buckling in the present system changes the Sis-Au-Sis 
angles to 101.75° (up-edge) and 101.29° (down-edge xl 
cell). Moreover, the angles between Au and Pb atoms 
are also slightly changed, and equal 81.78° (81.54°) in 
up-edge (down-edge) x 1 cells, to be compared with equi- 
librium value 81.97°. All this is reflected in the band 
structure, where very small energy gap develops in dis- 
persive Au-Pb-Si band (filled squares in Fig. [5|), when 
the step-edge buckling takes place. However, this small 
energy gap is slightly below the Fermi level. Similar ef- 
fect, although with slightly larger energy gap, has been 
found in Si(557)-Au surface [44], and has been assigned 
to the shortcomings of the LDA functionals, which fail 
to describe excitation spectra. So one can expect that 
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better description of the exchange and correlation effects 
(hke GW approximation) wih move the empty states to 
higher energies, thus driving the system into insulating 
state [45*, '46^. However, it is also possible, that present 
LDA results are correct, and the system stays in metallic 
phase. The step-edge buckling may not be able to open a 
true gap in Au-Pb-Si band, as it gives much smaller val- 
ues of Au-Si angles than in case of the Si(557)-Au surface. 
This last scenario seems to be supported by experimental 
data of Ref. [14], shown in Fig. [8] together with calcu- 
lated surface band structure for the step-edge buckled 
Si(335)-Au/Pb reconstruction. As one can see the agree- 




r K M' 



FIG. 8: (Color online) Comparison of the measured ARPES 
intensity (Ref. [14]) and the calculated surface band structure 
for the model with the buckling of the step edge. 

ment between measured and calculated band structure is 
very good. In particular there is no band pinned at the 
Fermi energy. 

To calculate the surface band structure shown in Fig. 
m following procedure has been applied. Since the bulk 
periodicity in direction [110] is 1 x a[i];o], and calculations 
have been performed in x2 unit cell in this direction, the 
resulting band structure features twice as many bands 
as in xl unit cell, folded back into first BZ of x2 struc- 
ture. This of course is not observed in experiment, as 
the ARPES is surface sensitive technique. To get rid of 
these bulk electron bands, only the surface atoms have 



been taken into account in calculations. The positions of 
these atoms have been obtained in full slab calculations, 
as described in Sec. [TTl To mimic rest of the Si layers, 
the surface Si atoms were passivated by hydrogen. The 
H atoms were placed in a way that the bond lengths and 
bond angles between surface and next from the surface Si 
layers were preserved. In this way calculated band struc- 
ture is shown in Fig. [51 Note that differences between 
in this way calculated surface band structure and that 
obtained from full slab calculations are negligible. 



VIII. CONCLUSIONS 

In conclusion, the structural and the electronic proper- 
ties of the monoatomic Pb chains on Si(335)-Au surface 
have been studied within the density functional theory. 
The obtained structural model features Pb chain on-top 
of the surface in the middle of terrace. As a result, two 
monoatomic chains are observed on single terrace, one 
made of lead, and the other one associated with the step 
edge Si atoms. Both chains are visible in STM experi- 
mental data. The STM topography is very well repro- 
duced within the present model. The calculated band 
structure shows clear one-dimensional character of the 
structure and remains in almost perfect agreement with 
the ARPES data. Finally, the buckling of the step edge 
Si atoms has been found, and unlike for the Si(557)-Au 
surface, does not drive the system into insulating phase. 
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